Character of eigenstates of the three- dimensional disordered Hamiltonian 
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We study numerically the character of electron eigenstates of the three dimensional disordered 
Anderson model. Analysis of the statistics of inverse participation ratio as well as numerical eval- 
uation of the electron-hole correlation function confirm that there are no localized states below the 
mobility edge, as well as no metallic state in the tail of the conductive band. We discuss also finite 
size effects observed in the analysis of all the discussed quantities. 
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I. INTRODUCTION 

Localization of electron in disordered systems,'^ man- 
ifests the wave character of the electron propagation. 
Components of the wave function, scattered on randomly 
distributed impurities, interfere with each other. This in- 
terference might lead to the electron localization. 

While all states are localized in one dimensional sys- 
tems, the localization of all electronic states in the 
three dimensional (3D) systems appears only when the 
strength of the disorder, W , exceeds certain critical value 
Wc- For weaker disorder, W < Wc, the energy band is 
divided into two parts, separated by the mobility edge 
Ec- It is supposed that all states with energy \E\ < Ec 
are metallic (conductive) , while only localized states exist 
for \E\ > Ec. 

Although electron localization is easy to understand 
intuitively, the wave character of electron propagation 
causes new, nontrivial phenomena in all three transport 
regimes. 

In the limit of weak disorder, the system is metal, 
but the scattering of electron on impurities is respon- 
sible for non-classical phenomena, such as the universal, 
system size independent conductance fluctuations^;- and 
weak localization (antilocalization)4i^ The complete de- 
scription of the transport was given by Green's function 
method,— random matrix theory,— DMPK equation^ or 
supersymmetric methodsi^ 

In the opposite limit of strong disorder, (localized 
regime), the fluctuations of the conductance are so 
strong, (they exceed the mean conductance in many or- 
ders of magnitude), that the conductance itself is not a 
relevant parameter of the theory any more. Instead, the 
logarithm of the conductance must be studied.— 

Owing to huge fluctuations, the analytical description 
of the critical regime is much more complicated. Thus, 
although the critical metal-insulator transition is well un- 
derstood by the single parameter scaling,— quantitative 
estimation of critical parameters is still almost unsolvable 
problem. Of particular interest is the critical exponent 
v^^ which determines the divergence of the correlation 
length, (X \E — Ec\^'^ &t the mobility edge. Over two 
decades, there is a discrepancy between theoretical pre- 



dictions and numerical estimations of v. For 3D Ander- 
son model, numerical results v « 1.5 — 1.57 obtained 
by various numerical method a^^d^d'^d^d^ , overestimate 
the analytical prediction, v = 1— li^. The disagreement 



is even worse in 4D (numerical data? 
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while theory predicts v — 0.5). Numerical data for v in 
low dimensional systems d — 2 + e (e <C l) ^^i^° also do 
not agree with the analytical e - expansions!^ 

In our opinion, the origin of this discrepancy lies in the 
procedure of the averaging over the disorderi2i2£ While 
averaging is extremely difhcult to perform analytically it 
is straightforward in numerical simulations. Therefore, 
we expect that numerically obtained critical parameters 
will be, soon or later, supported by analytical resultsj ^^i^^ 

The inhomogeneous spatial distribution of electrons in 
the critical and localized regimes^^'^^ inspire people to 
build mean field theories on the the analysis of the statis- 
tics of the local density of statej ^^'^^ Recent numerical 
data^^i^ led to new analytical theories of the transport 
in the strongly localized regime. —i^i^ 

In this paper we study numerically the character of 
eigenstates of the disordered 3D Anderson Hamiltonian. 
Our aim is to exclude any possibility for the existence 
of localized states below the mobility edge, Ec. The 
performed analysis is inspired by recent analytical the- 
ory of the Anderson transition^ which predicts that the 
number of metallic states decreases continuously to zero 
when Fermi energy approaches the mobility edge from 
the metallic side. The idea is formulated in terms of the 
electron- hole correlation function Tq{E,ijj) (defined later 
by Eqs. [3] and [3|) which possesses in the limit of small 
energy difference u and small wave vector q the diffusive 
pole of the form 



Tq{E,uj) = 



2ttp{E) 



'iA{E)uj + D{Lu)q'^ 



(1) 



{p{E) = Im G{E + ie)/TT is the density of states'^'' deter- 
mined by the one electron Green function G{E + ie), and 
D is the diffusion constant). 

Expression ([T]) differs from the "classical" diffusion 
pole^^ by the presence of the function A{E). It is claimed 
in Ref.'^'^ that A{E) increases when E approaches the mo- 
bility edge and becomes infinite at the critical point. The 
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ratio p{E) / A{E) determines the portion extended (diffu- 
sive) states from all available states, given by p{E). The 
rest states, p{E) x {A — 1)/^, are spatially localized, al- 
though E lies in the metallic phase, \E\ < Ec- 

Intuitively, the existence of localized states in the 
metallic phase seems to be impossible^. It also contra- 
dicts analytical analysis of the electron eigenstates^ i'^^'^^ . 
Nevertheless, no numerical analysis of this problem has 
been performed yet. The present paper fills this gap. 

We investigate in Section II the singular behavior of 
Tq^Ejiu) ~ uj~^ for q = and prove that A{E) = 1 for 
all energies E, both in the metallic and localized regime. 
This confirms theoretical expectationsiZi^. We also ana- 
lyze the diffusive pole in the metallic phase (band center) , 
find diffusive constant and discuss statistical properties 
of the function rq{E,uj). 

Another proof of the absence of localized states in the 
metallic phase is given in Section III, where we study the 
probability distribution of the inverse participation ratio 
(IPR)- defined later by Eq. [TT] Statistical properties of 
IPR were analytically studied in^ i35,36 ^ Statistical prop- 
erties of IPR at the critical point were subject of analyti- 
cal and numerical analysis in connection to the multifrac- 
tal spatial distribution of electrons Scaling of IPR 
in the critical region was proved imM. Here, we discuss 
how the probability distribution of IPR depends on the 
system size and the distance E — Ec oi the energy from 
the mobility edge. Our data show that the probability 
to find the localized state in the metallic phase decreases 
exponentially when the size of the system increases. 

Electron eigenenergies and wave functions are calcu- 
lated for the 3D Anderson Hamiltonian, 



(2) 



Here, r determines the site in the 3D lattice of the size 
L^, ep is the random energy distributed with the Gaus- 
sian distribution, Paiep) = {2ttW^)-^/^ exp{-el/2W^). 
Parameter W measures the strength of the disorder and 
V — 1 determines the energy scale. For E — 0, the 
critical disorder Wc ~ 6.15. We fix the strength of the 
disorder W = 2 throughout this paper. Then, the mobil- 
ity edge Ec — 6.58 separates the metallic and insulating 
phasei^ 



II. THE ELECTRON-HOLE CORRELATION 
FUNCTION 

In this section we investigate the electron-hole correla- 
tion function defined as 

r{E,Lj,r',r) = {G{E+Lj/2+ie,f,r')G{E-Lj/2-i€,?,r)). 

(3) 

Here, G{E + ie) = [E + ie - H] is the one-particle 
Green's function^, which determines the density of 
states, and (. . .) means averaging over realization of the 
disorder. 




o-o e = 5/L^ L = 18 
A-A e = 60/L^ L = 18 
□--o e = 80/L^ L = 18 
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FIG. 1: (color online) Probability distribution P(Im To) for 
E = 7, ui — 0.14, L = 18 and for three values of e. The density 
of states, p{E = 7) « 0.0076 and the mean level spacing is 
^ 130/L^ ^ 0.0223. The distribution obtained for e = 
is clearly unphysical, but the choice e — 60 /L^ is already 
sufficient to reach the Gaussian distribution of Im Fq. 



We calculate the Fourier transformation, 

T,{E,Lu) = ^e^'~(^~-'')r(i?,L.,P,r), (4) 

f,r' 

set q — and analyze the singular lj dependence 

B{E) 



ToiE.u;) 



(5) 



Comparison with Eq. ([!]) gives B{E) = 2ttp{E). Coeffi- 
cient B{E) would either equal to 2iTp{E) (case A = 1), 
or it would decrease to zero, B{E) ^ {Ec — E)", if the 
scenario proposed in Ref42. is true. 

Before presenting numerical data, it is worth to com- 
ment the numerical method of calculation of the Green's 
functions. Our method is based on the numerical inver- 
sion of the matrix E± ie — TL. For the reliability of data, 
the choice of the value of the small imaginary part of the 
energy, e, is crucial. We expect that e should be com- 
parable to the typical level spacing, l/{p{E)L^). Using 
numerical data for the density of state at the band cen- 
ter, p{E = 0) « 0.115, we choose e — This value 
is sufficiently large to avoid any numerical instabilities 
(discussed later in Sect. IIIC) in the band center, but it 
might be too small in the band tail, where the density of 
state is much smaller. Therefore, various larger values of 
e were used to guarantee the numerical stability of our 
results. 

As an example of how the value of e infiuences the 
accuracy of numerical results, we show in Fig. [T] the 
probability distribution P(Im Tq) for energy E = 7 and 
the system size L = 18. Statistical ensemble of = 2000 
samples was used to calculate the distribution. If e is 
much smaller than the level spacing, then the distribution 
P(Im To) consists of high peak close to zero, and very 
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FIG. 2; (Color online) The imaginary part of ro(i?,aj), given 
by Eq. ([5)l as a function of uj~^ . E = 5, L = 12 (circles), 
E — 6.58, L = 16 (squares) and E = 7, L = 18 (triangles). 
e — (open symbols). Higher values of e were used to 

check the stability of data: e = {E ^ 6), {E = 

6.58) and 80/L^ {E = 7) (full symbols). Only data left of the 
vertical line were used for calculation of B{E). Solid lines are 
fits In Im To = - In lj + In B(£). 
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FIG. 3: (Color online) The coefficient B{E) vs energy E com- 
pared with the density of states 2np{E). Only critical region 
is shown. Deep in the metallic phase, the coincidence of B{E) 
and 2np{E) is even better. Size of the system is L = 8 (circles) 
and L = 16 (triangles). Open triangles show B{E) calculated 
in Fig. [5] for e — 5/1/^ (triangle left) and for higher values 
of e (triangle right). The density of states is calculated for 
A'stat = 1000 systems of the size L = 16 (dashed line) and 
L — 8 (dotted line). Dirichlet boundary conditions were used 
for calculation of both p{E) and Fq. Vertical line marks the 
position of the critical point. 



Tq become reliable only for larger values of e, for which 
the density of states is a smooth function of the energy. 
As shown in Fig. [U P converges to the Gaussian distri- 
bution, independent on e for sufficiently large values of 



A. Singularity of Fq for a; — > 

Figure [2] shows numerical data for the imaginary part 
of Fo(-E,ijj) as a function of u; for three values of the 
energy E — 5, E ~ 6.58 (the mobility edge) and E = 7. 
Data prove the singular behavior Fq ~ 1/uj. Comparison 
of numerical data calculated for two and more different 
values of e enable us also to estimate the accuracy of 
our results. Although the singularity ^ l/w transforms 
to Im Fq ~ + e^) when e ^ 0. Im Fq becomes 

independent on e for u; ^ e (the region left from the 
vertical line in Fig. [J). Only these data were used for 
the calculation of the coefficient B{E). 

In Fig. [3] we plot numerical data for the coefficient 
B{E) and compare them with 2ttp{E). The density of 
state p{E) was calculated by direct diagonalizing of the 
Hamiltonian for L = 8 and L = 16. To increase the 
number of eigenstates, statistical ensembles of = 10"^ 
samples were used for each L. As shown in Fig. [31 the 
density of states in the band tail still depends on the 
system size. Clearly, L = 8 is not sufhcicnt for the cal- 
culation of p. To check the convergence of the density of 
states, we calculated the density of states for the energy 
E = 7, also from the statistical ensemble of samples of 
the size L = 40. The eigenenergies were calculated by 
Lanczos algorithm. A comparison of the obtained den- 
sity of states for L = 16 and 40 [p{E = 7) = 0.0076 for 
L = 16, and p{E = 7) = 0.0087 for L = 40] indicates 
that the convergence of the density of states is very slow 
in the band tail. 

As shown in Fig. [31 obtained coefficient B{E) agrees 
very well with our data for the density of states, 



B{E) = 2ttp{E). 



(7) 



The agreement is even better when we compare p and 
B{E) calculated for the same system size. Since the size 
corrections of both B{E) and p{E) are positive (both 
quantities increases when L increases), we conclude that 
our data for B{E) showed in Fig. [31 definitely exclude 
the possibility that B{E) decreases to zero when E ap- 
proaches the mobility edge. 



long tail towards high values. This is because the density 
of states 

consists of set of very narrow separated peaks centered 
around eigenenergies En for small e. Numerical data for 



B. Fg for non-zero q 

In this section we calculate the electron-hole correla- 
tion function Tq for non-zero values of q. We show that 
numerical data are consistent with theoretical prediction. 
From numerical data, we estimate the diffusion constant 
D. 
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FIG. 4: The real and the imaginary part of q'^rq{E,Lu), given 
by Eq. as a function of u/q^ for E = Q (the band center) 
and for three values of <j, q = (27r/L)[l, 0, 0], (27r/L)[l, 1, 0] 
and (27r/L)[l, 1, 1]. Solid lines show theoretical prediction, 
given by Eq. (|8]). with D = 1.055. Dashed lines are fits of 
numerical data for L = 16 with 0.94 < D < 1.15. 
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FIG. 5: (color online) The probability distribution of q^Tq for 
q — 1 and three system size, L — 12, 14 and 16 and Lu/q'^ ~ 1. 
The width of the distribution does not depend on the size of 
the system. The mean value, (q'^Vq) = 0.37, and the variance, 
var q^Tq « 0.11. does not depend on the size of the system. 



and (1, 1, 1) (in units of 2tt/L). As discussed above, e — 
is already sufficient for numerical analysis of the 
electron-hole correlation function at the band center. We 
use this value in all calculations below. 

Numerical data lie on the one universal curve. This 
universality is better pronounced for small values of q. 
Stronger finite size effects are observed for larger q. Data 
for L = 16 are fitted to Eq. ^ with fixed density of 
states, 2ttp{E = 0) = 0.72 and free parameter D. Fits are 
shown in Fig. lUby dashed lines. From fits, we estimate 
the diffusion constant. 



D « 1.05 ±0.10. 



(9) 



This value is compared with diffusion constant calculated 
by the transfer matrix method from the dependence 
of the conductance, 



In general, D is a function of both uj and q. Numeri- 
cal analysis of D{uj, q) for critical disorder W — Wc and 
energy E = Q was performed in Ref."^^ Numerical simula- 
tions confirmed scaling behavior of the diffusive constant 
in the critical regime, predicted theoretically Since the 
critical region around the mobility edge E ^ Ec is very 
narrow,— we restrict our analysis to the metallic regime, 
W = 2 <^ Wc and E — 0. Here, we expect that D is 
constant, independent on the frequency and wave vector. 

Fig. [3] shows the real and the imaginary parts of the 
function q^Tq, 



g(L) = aL 



(10) 



2t^ 

1 



{E^u) 



2ttp{E) 
-iuj/q^ + D' 



(8) 



as a function of w/g^. The size of the system is L — 12, 14 
and 16 with periodic boundary conditions. Three values 
of the wave vector were considered: q: (1,0,0), (1,1,0) 



where a = e^D{E)p{E) is the conductivity's We ob- 
tained D « 1.055, which perfectly agrees with our esti- 
mation 

Finally, we present in Fig. [5] the probability distri- 
bution of real part of q^Tq for three sizes of the system 
bj/q'^ Ki 1 for all systems. Our results confirm that the 
distribution depends only on the ratio lojq^. This is con- 
sistent with Eq. ([5]). More important, data in Fig. O 
indicate that P does not depend on the system size. 



III. INVERSE PARTICIPATION RATIO 

The absence of localized states in the metallic phase 
can be demonstrated also by the analysis of the proba- 
bility distribution of inverse participation ratio, defined 
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/(i?„) = ^|$„(r)| 



(11) 



Here, En and $„(r) is the nth eigenenergy and eigenfunc- 
tion of the Hamiltonian respectively. If the \E\ < E^, 
then we expect that aU eigenstates are conductive, with 
the wave functions distributed moreless homogeneously 
throughout the sample, so that |<i>„(r)p oc L^"*. In- 
serting in Eq. (fTT|) we obtain that /(£"„) oc L^^{m 
the 3D system). On the other hand, the wave function 
of localized electron is non-zero only in a small region, 
where |$n(T')| ^ 1- Hence, I{En) 1. The size de- 
pendence of I (En) in the critical region deserves more 
detailed analysis since the spatial distribution of electron 
is multifractal^^'^"^ and / oc L^'^^ where d2 ~ 1.28^ 



A. The size dependence of IPR 

The energy spectrum of the Hamiltonian depends on 
the system size, L, and on the microscopic details of the 
disorder in a given sample. For a given system size, we 
consider a statistical ensemble of Ns samples which dif- 
fer only in the realization of random energies, ep. For 
each sample, we calculate all eigenenergies, i?„, lying in 
a narrow energy interval, E — 6, E + S, and calculate cor- 
responding / = I (En). For the ith sample, the number 
of eigenstates, rii, depends on the microscopic realization 
of the disorder. 

Collecting A'stat — J2i values of IPR, we can con- 
struct its probability distribution P{I) or P(ln/). Since 
the values of / might fluctuate in many orders in magni- 
tude in the critical region^, it is more convenient to use 
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FIG. 6: (Color online) The left upper panel shows the sys- 
tem size dependence of y oc a_ElnL with as = —3.006, 
aa = —2.326 and ae.s = —1.48. Three other panels show 
the probability distribution P(ln/) for three energies of the 
electron and for the system size L = 18 (v)i 24 (A), 32 
(<>) and L = 40 (o). Note the scaling of the horizontal axis. 
A^stat = 10^ eigenstates from the interval \En — E\ < 0.025 
were used for the statisticsi^ 
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FIG. 7: (Color online) Left panel: The size dependence of 
the probability H, defined by Eq. with 7min = — 21nL 

for the energy E — 3. The right panel shows the probability 
n^iiin lor -^min = —3/2 In L and E = 6. Data confirm that 
Hj^min decreases exponentially when the size of the system 
increases. 



the logarithm of / and the mean value. 



Y = Y{E) = 



■Y, J2 ln/(i?«). (12) 

i |E-£;„|<<5 



The upper left panel of Fig. [S] shows the system size 
dependence of Y for three energies below the mobility 
edge. We found that y ~ a^; InL. In the metallic regime, 
we expect = —3 for all < ii^c- At the mobility 
edge, aE=E^ — ^2, the fractal dimension^^. For E — 
3, which is the energy far below the mobility edge, we 
indeed find ay, — —3, in agreement with our expectations. 
Higher values of a^; obtained for energies closer to the 
mobility edge are due to the finite size effects. We expect 
that aE converges to —3 when the size of the system 
increases, L — > oo. This is consistent with analytical 
expression for the mean IPR^ 



{I{L))=L-^[l + ^L/m. 



(13) 



Deep in the metallic regime, \E\ < Ec, the conductance g 
is oc L (Eq. Iini), so that Eq. ^ reproduces (/) ~ 
However, the linear increase oi g L can be obtained 
only when the size of the system L 3> ^. For smaller-size 
system, L the correction term 4:L/£g in Eq. (fT5)) 

becomes L dependent and causes the deviation from the 
L~^- dependence of mean /. The scaling behavior of IPR 
for energies close to the mobility edge is discussed later 
in section UlI CI 



B. Probability distribution of IPR 

Because of the randomness and wave character of the 
electron motion, the mean value of any quantity might 
not provide us with the entire information about the sys- 
tem. For instance, although Y oc — 31nL, there still 
might exists some localized electronic states with an 
eigenvalue £"„ and In/ ~ 0. To measure the proba- 
bility that insulating states exist in the metallic phase 
{\E\ < Ec), we plot in Fig. [Hlthe probability distribution 
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FIG. 8: (Color online) Probability distribution P(ln/) for the 
localized states for L = 12, 16, 20 and 24. The distribution 
only weakly depends on the size of the system and decreases 
exponentially when In / decreases. Consequently, the proba- 
bility to find, for instance, the eigenstate with In/ < — InL 
decreases exponentially when L increases. Therefore, the 
probability to find the conductive state in the localized phase 
(|/5| > Ec) is exponentially small. 



P(ln/) calculated for three energies and various size of 
the system. 

Data for E = i confirm that P(ln/) gets narrower 
when L increases. This is consistent with analytical re- 
sult, var / ~ L~^^ The narrowing of the probability 
distribution is less visible for energies close to the mobil- 
ity edge Ec- 

Since the mean value, 1", decreases as ^ aE^T^L when 
L increases, the existence of localized states is possible 
only if P(ln/) possesses a long tail which assures a non- 
zero probability to have In/ for any system size. 
However, our data in Fig. [Slshow that this is not the case. 
Contrary, P(ln/) decreases exponentially for In / larger 
than Y . To measure this exponential decrease quantita- 
tively, we calculate the probability that In / is larger than 
certain value, Imm- 



dln/'P(ln/') 



d/'P(/') (14) 



for ^; = 3 and /„ 



L-3/2 for 



We choose /,„] 
E = 6.0. In Fig. [7] we prove that Hi^^^ decreases expo- 
nentially as a function of the size of the system L. This is 
consistent with theoretical prediction P(/) ~ exp(— a/)j^ 
Note, this exponential decrease is visible already for en- 
ergies very close to the critical point, (right panel of Fig. 
[7]). Since /(/?„) ~ 1 for the localized state the prob- 
ability to observe the localized state inside the metallic 
phase decreases exponentially when the size of the sys- 
tem increases. We conclude that the probability to find 
any localized state is zero in the limit of L — > cx). 

Similar conclusion, namely that there are no metal- 
lic states in the energy interval E > E^ can be derived 
for localized phase. In Fig. [5] we show the probabil- 
ity distribution P(ln/) for eigenstates around the energy 
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FIG. 9; Probability distribution P(ln/) for various energies, 
and three sizes of the system: L = 12 (top), L — 16 (middle) 
and L = 24 (bottom). Shaded area is the critical distribution 
Pc for E — Ec- Note, horizontal axis is scaled by —da In L 
in order to keep the critical distribution for Ec = 6.58 in the 
center of the figure. 



E = 7.5. Clearly, the distribution is size independent, 
and decreases exponentially for small values of In /. Since 
metallic state requires that In/ ~ —3 InL, we conclude 
that there are no metallic states in the insulating phase. 



C. Finite size scaling 

The right upper panel of Fig. [6] shows the distribution 
of IPR in the metallic regime. The distribution is cen- 
tered at 3 InL, in agreement with our expectation. More 
important is the form of the distribution for larger val- 
ues of In /. Our data show that the probability to observe 
In/ —2 InL decreases exponentially when the size of 
the system increases. Hence, we conclude that there are 
no localized states. 

However, as is shown in lower panels of Fig. [51 nar- 
rowing of the distribution P(ln/) can be numerically ob- 
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served only when the energy E hes deep in the metalhc 
phase. This is easy to understood. The metalhc phase 
can be observed only in systems of size i :» C {E) , where 
^{E) is the correlation length. Since £,{E) diverges as 
^{E) cx \Ec — E]^" iv w 1.57 is the critical exponent) 
when E — > E~ , we cannot observe the metallic behavior 
for energies close to the mobility edge and for fixed size of 
the system. Nevertheless, reliable conclusion about the 
character of metallic states in the vicinity of the mobility 
edge can be drawn with the use of the finite size scaling 
analysis>i^ The probability distribution P(ln/) depends 
both on the energy E and on the size of the system L. 
Following the scaling theory, P calculated for a given en- 
ergy E and size L is equivalent to P obtained for E' 
closer to Ec, but larger system size L' > L. 

Scaling of IPR in the critical region was numerically 
proved in Ref.'^'*. In Fig. [9] we demonstrate how the 
scaling idea works in the metallic phase. The distribu- 
tions P(ln/) were calculated for various energies of the 
electron and for three size of the system. We see the simi- 
larity in the form of P(ln J) calculated for different E and 
L. For instance, the distribution for E — 3.1 and L = 16 
is similar to the distribution for E — 0.1 and L — 12. 
Similarly, we can compare i? = 5.1 and L = 24 with that 
for E = 4.1, L ~ 12. We also see that the form of the 
probability distribution P(ln/) only weakly depends on 
the energy E of the electron when L ^ ^{E). From this 
similarity we conclude that the properties of the distri- 
bution P(ln/) are universal in the metallic phase when 
L — > oo. Therefore, the probability to find the localized 
state decreases to zero for any energy \E\ < Ec- 

IV. CONCLUSION 

We studied numerically two parameters, important for 
the construction of the analytical theory of the metal- 
insulator transition. First, we verified the relation be- 
tween electron-hole correlation function F and the den- 
sity of states, Eq. ([5|). We proved that this relation not 
only holds for all energies of the electron, both in the 



metallic and localized phase, but can be recovered for 
any size of the system. Our numerical procedure enables 
us to calculate, from F, the diffusion constant D. In the 
metallic limit, D agrees with estimation from the trans- 
fer matrix. Also, numerical data indicate that Tq{E,uj) 
is not the self-averaged quantity in the metallic regime. 

It would be interesting to investigate also the scaling 
behavior of the diffusive constant in the critical regime4S 
Such analysis could confirm numerical scaling observed 
recently for the case of the critical disorder at band 
center. However, since only frequencies to > e are rele- 
vant in numerical data, we have to analyze much larger 
system size in order to fit both energies E ±uj/2 into the 
narrow critical region. ■^^ 

We also present numerical data for the mean values 
and for probability distributions of the inverse participa- 
tion ratio. Our data, consistent with previous analytical 
results, enable us to prove that there are no localized 
states inside the metallic phase. All electron states are 
extended, and the probability to find the state which does 
not span through the sample decreases exponentially to 
zero when the size of the system increases. Similarly, 
no metallic states were observed on the opposite side of 
the critical point: in the insulating tail of the spectra, all 
electronic states are localized. Although this result seems 
to be easily accepted^, it was never proved numerically. 

In contrast to the analytical theory, numerical methods 
do not enable the analysis of the behavior of any quantity 
in the limit of infinite system size. We can only describe 
how the variables of interest change when the size of the 
system increases. With the use of the finite size scaling 
hypothesis, we conclude that our results remain valid also 
in the limit of L — > oo. Since all data were obtained 
without any additional assumption about the averaging 
procedure or the statistics of parameters of interest, they 
can serve as a starting point for the construction of the 
analytical theory of the Anderson transition. 
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